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Abstract 

Local porosity distributions for a three-dimensional porous medium 
and local porosity distributions for a two-dimensional plane-section 
through the medium are generally different. However, for homoge- 
neous and isotropic media having finite correlation- lengths, a good 
degree of correspondence between the two sets of local porosity distri- 
butions can be obtained by rescaling lengths, and the mapping associ- 
ating corresponding distributions can be found from two-dimensional 
observations alone. The agreement between associated distributions 
is good as long as the linear extent of the measurement cells involved 
is somewhat larger than the correlation length, and it improves as the 
linear extent increases. A simple application of the central limit the- 
orem shows that there must be a correspondence in the limit of very 
large measurement cells, because the distributions from both sets ap- 
proach normal distributions. A normal distribution has two indepen- 
dent parameters: the mean and the variance. If the sample is large 
enough, LPDs from both sets will have the same mean. Therefore 
corresponding distributions are found by matching variances of two- 
and three-dimensional local porosity distributions. The variance can 
be independently determined from correlation functions. Equating 
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variances leads to a scaling relation for lengths in this limit. Three 
particular systems are examined in order to show that this scaling 
behavior persists at smaller length-scales. 
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1 Introduction 



The study of porous and heterogeneous media is scattered throughout many- 
fields of science, including mathematics, solid state physics and chemistry, 
materials science, geology, hydrology, environmental technology, petroleum 
engineering and separation technology. The central problem is the specifica- 
tion of the random microstructure, which is needed to predict macroscopic 
physical properties ( see |]ll0[000 for overviews). 

A complete specification of the random microstructure is both impracti- 
cal and unnecessary. It is therefore important to have a general statistical 
description of the microstructure available. Such a description should meet 
four criteria: it should be well-defined in terms of geometrical quantities, 
it should involve only experimentally accessible parameters, it should be of 
economical size and it should be usable in exact or approximate solutions 
of the underlying equations of motion. Currently there are only two sta- 
tistical methodologies available which fulfill all four requirements; these are 
correlation functions [| § g] [| fig in| [|l3l H |5[ and local geometry dis- 
tributions [0 ago 00. 

The recently developed local geometry distributions are a functional gen- 
eralization of the correlation-function approach]^. The aim of the present 
paper is to increase the practical applicability of local geometry distributions. 
Interest in these distributions arises from their potential use in distinguishing 



between different microstructures[|T^]. Increasing their practical applicability 
is particularly important because local geometry distributions can be used 
in mean field theory calculations of the dielectric and transport properties of 
various porous media, for instance porous rock filled with brine [|l^ ITSf. Pre- 



liminary comparisons of the theory with experimental results have already 
been carried out, and show good agreement p!9| pO|| . 

Two main objectives will be pursued in the present paper. First, we 
show the feasibility of obtaining three-dimensional local porosity distribu- 
tions(LPDs), by presenting data from a real sandstone specimen. Second, we 
show to what extent three-dimensional LPDs can be obtained from purely 
two-dimensional observations. The dependence of LPDs upon the dimension- 
ality of the system is a consequence of their dependence on a measurement 
cell or observation region over which geometric observables, such as porosity, 
are averaged. The possibility of extracting three-dimensional distributions 
from two-dimensional data arises from limit theorems which guarantee that 
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the asymptotic LPDs become independent of the shape of the measurement 



We begin our discussion with a definition of LPDs. Next we present a 
simple argument based on the central limit theorem showing that for large 
measurement cells, three-dimensional LPDs can be mapped onto two-dimen- 
sional LPDs by a simple rescaling of lengths. Then, we show that this map- 
ping becomes exact for the percolation model. Subsequently, we analyse the 
overlapping sphere (or continuum percolation) model, where the mapping 
holds only asymptotically. Finally we present the first measurement of three- 
dimensional LPDs. These are obtained from the pore-space reconstruction 
of a real sandstone specimen. 

The main result of our investigation is the fact that the mapping between 
two- and three-dimensional LPDs can be determined from purely two-dimen- 
sional observations if the medium is homogeneous and isotropic. For general 
anisotropic media however, the family of three-dimensional LPDs can only 
be determined from a three-dimensional measurement. 

2 Basic Concepts and an Expression for Vari- 
ance 

The discretized porous medium is considered to occupy a set S C 7/ ^ where Z 
is the set of integers. Pore-space is denoted by P, where P C S. Considering 
for simplicity only two-component media, we define matrix-space M as the 
complement of P in S, M = S — P. The characteristic function of a set A is 



celli. 




1, rGA 
0, otherwise 



where r e Z''. To describe the state of each point in the system, we label 
each point in S with an index i and define the random variables 



X, = X(ri)=xp(ri) 



A specific geometry is then specified by the A^-tuple 



^ = (Xi,...,X;v), iV = |S 
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where |S| is the number of elements in S. 

The probabihty space for the system is the triple {G,S,P). G is the 
sample space consisting of all possible geometries g. The event space £^ is a 
Boolean algebra of subsets of G. The probability P of finding a geometryQ^f is 
a finitely additive non- negative set function defined on S, such that P(0) = 
and -P(G) = 1. A stochastic medium is homogenous if P is invariant under 
translations. Isotropy is defined analogously as invariance under rotations. 
The expectation value of a random variable f{g) is 

(fig)) = E fi9)Pi9) 

geG 

and in the following, we will use the notation 

n 

^^(ri, . . . ,r„) = (nXp(ri)) 

i=l 
n 

Cn{ri, . . . ,Yn) = (nxp(r*) - (Xp(ri))) 

i=l 

for the moment functions S'n(ri, . . . , r^) and cumulant functions C„(ri, . . . , r^) 
of order n. For homogeneous media, 5*1 (r) = {(f)) is the expected bulk poros- 
ity. 

Given a measurement cell denoted as K (K C S), we define the local 
porosity of K as the random variable 

II II j.gz'' 

The local porosity distribution(LPD) [p^ for K can be defined as0 

where 5(0 — 0(K)) is the Dirac delta function, and the first moment of this 
LPD is 

Jo 

^Strictly speaking, P is defined on sets of geometries E e f . However, we assume tfiat 
{g} e £ for all geG. 
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Using the definition of /i, and clianging tlie order of summation gives 



geG 



E Xp(r)XK(r)\ = p^E(^pW)^K(r) 



wliere we liave also used the definition of 0(K). Thus, we obtain the expres- 
sion 

J(K)= /VM0;K)c?0=(0(K)) = -^5:Si(r)xK(r) (1) 



For homogeneous media, Si = (0), the bulk porosity, so that (K) = 
(0). Similarly, we can show that the local porosity variance far[0(K)] = 

lo (0-0 (K)) /i(0; K) d(j) is given by 

var[m)] = T^ E C2(r,r')xK(r)xK(r') (2) 



r.r'eZ 



d 



If the medium is homogeneous we can perform one of the sums. In this 
case C2(r, r ) = C2(0,r — r). Defining y = r — r, and using the fact that 
XK(r + y) = XK-y(i"), where K — y is the set K translated by the vector 
— y, we obtain 

var [0(K)] = '^^|1^2~^^' g'2(0> y) (3) 

Note that the variance var[(f)(K)] is simply related to the granularity con- 
cept which is used for photographic materials [pi]] . The relations above are 
applicable to infinitely large realizations of homogeneous stochastic media. 



3 The Central Limit Theorem and Scahng 

Consider again the class of homogeneous media in which spatial correlations 
of the random geometry decay with a finite correlation-length ^. For such sys- 
tems, we can apply the central limit theorem to measurements of local porsity 
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for very large measurement cells. To do this, we partition a rf-dimensional 
hypercubic measurement cell K of side L into n = [L/aY hypercubic subcells 
{Kj}, i = 1, 2, . . . , n, of side a so that 



K= IJK^, K,nK,=0 for i^j 



and rewrite 0(K) as the average 

over all the subcell porosities 0(Kj). If we choose L and a such that ^ ^ a <C 
L, then n 3> 1, and the central limit theorem becomes applicable, because 
the 0(Kj) are weakly correlated random variables. The central limit theorem 
ensures the existence of a number a such that 



M0;K)^t^^exp 
v27r(T 



(0 - 4>f {L/aY 

2(72 

were the accuracy of the approximate equality improves as, L/a increases. In 
the limit L/a ^ oo subject to a ^ ^, one has ^ 6{(f) — 0), independent 
of dimension and independent of K. Hence, for L/a sufficiently large but 
still finite, the form of ii{(f>; K) depends only on the bulk porosity (p and the 
variance t'ar[0(K)] = a^/ [L/aY, t)ut not on the shape of K. If the medium is 
also isotropic, then a sufficient condition for the equality of two such limiting 
distributions for d and d -dimensional measurement cells in the same medium 
is thus 

(L'\ (a'V^'^ (L\'l^ 



o j \a) 

because will be the same in both cases. 

Thus we find that for homogeneous isotropic systems with finite correlation- 
length, two- and three-dimensional LPDs for large L are simply related by a 
rescaling of lengths. We now examine the site percolation model, for which 
the scaling relation is exact. 
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4 The Percolation Model 



In the site percolation model, each lattice site has a probability (0) of being 
occupied and (1 — {</))) of being empty that is independent of the other 
sites. Therefore, if we consider occupied sites to be pore and empty sites as 
matrix, /i(0;K) depends only on (0) and |K|, and is given by the binomial 
distribution. If the underlying lattice is hypercubic, then the LPD for a 
d-dimensional hypercubic measurement cell K of side L is 

where 0j = i/L'^ and i = 0, 1, ... , L'^. Measurement cells consisting of the 
same number of points must have identical LPDs. Therefore, if K is a 
d -dimensional measurement cell of side L 

K') =/x(0,;K) if L'"^' =L'^ 

The scahng relation 

L' = L'^/'^' (5) 

is not solvable in Z for all L, L' . Nevertheless, there exist infinitely many 
solutions, since each positive integer i yields the solution L' = i'^ and L = i'^ 
Therefore, for the site percolation model, two- and three-dimendional LPDs 
are, within the constraints imposed by the discrete nature of the model, 
simply related by a rescaling of lengths. 



5 The Overlapping Sphere Model 

An often useful grain model for porous media is obtained by randomly dis- 
tributing overlapping spheres. The spheres represent grains of matrix ma- 
terial. A model configuration is typically created by randomly distributing 
point-centers with constant point-density in a given volume and attaching 
a sphere of radius ro to each point-center. The system analyzed below 
was obtained by randomly distributing 8000 point-centers in a volume of 
256x256x256 voxels. A value of tq = 9.45 voxels has been used, and periodic 
boundary conditions have been implemented. The bulk porosity of the over- 
lapping sphere model is (0) = exp —p^nrQ , where p is the point-density. 
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The point density and sphere radius were chosen such that (0) = 0.1855, 
matching the porosity of a real sandstone specimen discussed below. A plane- 
section showing the degree of discretization is displayed in Figure |I[ 

Three-dimensional LPDs for L G {1, 2, . . . , 64} have been measured, and 
the resulting family of distributions is shown in Figures ^ Two-dimensional 
LPDs were obtained for L G {1, 2, . . . , 150} by averaging the measured dis- 
tributions for all cross-sections perpendicular to the z-axis. The results are 
displayed in Figure ^ To investigate the degree of correspondence between 
two- and three-dimensional LPDs we calculate the quantity 

/■I 

A(L, L) = I /i3d(0; L) - fi2d{4>] L)\d(f) 
Jo 

as a quantitative measure of the deviation between two distributions. For 
fixed L, we minimize A(L, L ) by varying L . The matched distributions 
selected by this criterion are shown in Figure ^ for a choice of (L, L )-values 
indicated in the figure for each matching pair. The agreement between the 
matched two- and three-dimensional distributions is satisfactory. The agree- 
ment for large L stipulated by the central limit theorem improves with better 
system-statistics. This is seen from the matched LPDs shown in Figure ^, for 
a system of 64000 instead of 8000 spheres, having the same voxel- dimensions 
and bulk porosity. The poorer sphere-resolution of ro = 4.73 pixels makes 
this configuration less suitable for analyzing /i at smaller length-scales, so we 
shall not discuss this configuration further. 

The matching of distributions based on A(L, L ) is shown as the solid line 
in Figure ^ Another matching of distributions is obtained from the criterion 
of equal variances established in section 3. The mapping L (L) obtained from 
the equal-variance criterion is shown as the dashed curve in Figure For 
lengths that are not too large relative to the size of the system, both criteria 
give more or less the same rescaling of lengths L (L). Note, however, that 
the discretization can affect L (L) at small L. 

In practical applications, it is important to estimate the three-dimensional 
LPDs from observations on two-dimensional cross-sections alone. In the 
matching above, however, it was assumed that the three-dimensional LPDs 
are known. A way out of this dilemma is provided by equations @ and (Q). 
The agreement between the dashed and the solid lines in Figure |^ validates 
the matching criterion based solely on the variances expressed in equation 
(^. On the other hand, equation allows us to estimate the variance from 



9 




Figure 1: A digitized plane-section from a 256x256x256- voxel digitized con- 
figuration of overlapping spheres, consisting of 8000 randomly distributed 
spheres, all of radius Tq = 9.45 voxels. 
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Figure 2: The family of three-dimensional LPDs ^3d{4'', L) for L G 
{1,2,..., 64} for the overlapping sphere model. Note that the position of 
the saddle-point is near L = 18.9, which is the correlation length for the 
system. 
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Figure 3: The family of two-dimensional LPDs pL2d{4>'-, L') for L' G 
{1, 2, . . . , 150}. Although this surface is not the same as the in the pre- 
vious one, the previous surface can be approximately recreated from this one 
by rescaling the L axis. 
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'(24, 39) 



Symbols: 




|J3d(*f 


;L) 


|J2d(*f 


;L') 




0.0 



0.2 0.4 0.6 

local porosity (j) 



0.8 



1.0 



Figure 4: A selection of matching two- and three-dimensional LPDs from 
the mapping obtained by minimizing A(L, L ). Solid lines indicate three- 
dimensional LPDs and dashed lines show two-dimensional LPDs. The side- 
length of the measurement cells for matched LPDs are shown on the figure 
as (L, L )-pairs. 
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0.0 0.2 0.4 0.6 0.8 1.0 

local porosity 



Figure 5: Selected A(L, L')-matched distributions, similar to those shown 
in the lower half of the previous figure, from an ensemble of 64000 spheres 
of tq = 4.73 voxels, within a volume of 256x256x256 voxels. The matching 
(L, L )-pairs are indicated in the figure. 
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Figure 6: The scaling relation L' (L) associating two- and three-dimensional 
LPDs, as determined by four different criteria. The solid line is based on 
minimizing A(L, L ). The remaining curves are based on equating variances: 
the dashed curve is obtained from the actual distribution variances, the dot- 
ted curve and the dash-dotted curves are obtained from variances computed 
from (H) using measured and analytical values for C2 respectively. 
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6*2 (0, r) measured on two-dimensional cross-sections as long as the system is 
isotropic. The results of such a matching procedure are shown as the dotted 
and the dash-dotted curves in Figure ^. The dotted curve uses a direct mea- 
surement of C2 from two-dimensional sections, while the dash-dotted curve 
uses the exactly known analytical form of C2. The direct measurement of 
C2 uses the same series of planar-sections as the measurement of the two- 
dimensional LPDs; the correlation function is first calculated in each plane 
and subsequently averaged over all planes. The necessity to average over 
all planes, and the influence of discretization effects were recently discussed 
systematically [^]. The resulting averaged two-point correlation function 



G(0,r) 



C2(0,r) S2(0,r)-Si(0)5i(r) 



C2(0,0) S2(0,0)-5i(0)5i(0) 



is shown in Figure ^ Because the model of overlapping spheres is homo- 
geneous and isotropic G(0,r) = G{r). The exact analytical expression for 
5*2 (0, r) is given as 



(0,r) = S2(|r|) = exp 



p is the density of point-centers, d is the sphere-diameter, V3(A) is the volume 
of the set A and B(r, c) = {x : |x — r| <c} is the ball of radius c around r. 
The exact G{r) vanishes for r > d, and is shown as the dashed curve in 
Figure ^. Visible discrepencies between the measured and analytical curves 
are due to insufficient system-statistics. 

The satisfactory agreement of the dotted and the dash-dotted curves in 
Figure ^ results from the isotropy of the overlapping sphere model. The 
isotropy condition G{r) = G{r) guarantees knowledge of the full function 
G{r) from a measurement along any ray. Knowledge of the full function is a 
prerequisite for calculating the three-dimensional variance from (|^). Hence, 
the above precedure is not applicable for general anisotropic media. 

To assess the relative reliability of our rescaling procedure for isotropic 
systems, we compare in Figure | the two- and three-dimensional variances 
computed from the exact and measured correlation functions G{r) with those 
computed from the local porosity distribution /i(0;K). Solid lines are vari- 
ances computed from yu(0; K), dashed lines are variances computed from the 
measured G{r) and dotted lines are variances computed using the analytical 
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10 20 30 

r 



Figure 7: The 2-point correlation function G as a function of distance r. The 
sohd curve shows measured values, while the dashed curve is the analytical 
result for an infinite sample. 
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G{r). The overall agreement is satisfactory, given that the side length of the 
whole system is only 256 voxels. 



6 A Sandstone specimen 

In this section, we analyse experimental data. The three-dimensional pore- 
space for a specimen of Savonnier oolithic sandstone was reconstructed from 
serial thin sections||2^. The data-set consists of 99 binary images, each of 
1904x1904 pixels, representing planar cross-sections through the medium. 
The distance between the planes of consecutive images is 10/im. This is 
also the pixel size. Each section was obtained by thresholding a digitized 
micrograph into a black and white image. The threshold value was chosen 
to match the experimentally-known bulk porosity of = 0.1845. Figure ^ 
shows part of a typical image measuring approximately 1 cm along each side. 



In Figure |I0| we show a three-dimensional cubic subset with a side length of 
1 mm, in which the pore-space is shown in black, the matrix is transparent 
and the pore-space boundary is gray. 



Figure |T0| shows that the pore-space is anisotropic. This is corroborated 
in Figure |lT], where the correlation function G\\{r), measured in the plane 
of the thin sections is shown together with G±{r) measured in the direction 
normal to the thin sections. The anisotropy is at least partly caused by the 
misalignments of individual images and variations in lighting and threshold- 
ing during the reconstruction. The correlation function G\\{r) also shows 
the existence of large-scale heterogeneities. Given these imperfections we 
do not expect to be able to reconstruct the three-dimensional LPDs from 
two-dimensional information alone. 



The family of three-dimensional LPDs is displayed in Figure |T2|. The fam- 
ily of two-dimensional LPDs is displayed in Figure [l^. The three-dimensional 
distributions have been measured on a sub-system of 1881x1881x99 voxels 
for measurement cell sizes L e {1, 2, . . . , 96}. The bulk porosity of the sub- 
system is 0.1861. The results displayed in Figure O are, to our knowledge. 



the first measurements of three-dimensional LPDs on a natural rock sample. 
The two-dimensional LPDs shown for L G {1, 2, . . . , 240} were obtained, as 
before, by averaging the distributions for each plane. 

Figure |l^ shows the two- and three-dimensional variances computed di- 
rectly from the distributions as solid lines. The dashed lines are variances ob- 
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Figure 8: Three different determinations of variance for botli two- and tfiree- 
dimensional LPDs. Tfie actual distribution-variances are tlie solid curves. 
The dashed and dash-dotted curves are computed from (H) using measured 
and analytical values for C2 repectively. 
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Figure 9: Part of a typical planar-section through the Savonnier oolithic 
sandstone. The side-length of the image is ~ 1 cm. The pore-space is shown 
in black and rock matrix in white. 
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Figure 10: A roughly cubic sub-volume of 100x100x99 voxels from the same 
sandstone shown in the previous figure. The pore-space is shown in black, 
matrix is transparent and the pore-space boundary is gray. 
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Figure 11: Two-point correlation functions for the sandstone specimen. 
G\\{r) and G±{r) are shown as the sohd and dashed curves respectively. 
These functions are shown as functions of r + 1 in order to plot the r-axis in 
logarithmic scale. Note the presence of several correlation-lengths in Gy . 
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Figure 12: The family of three-dimensional LPDs fi3d{4>] L) for L e 
{1,2, ...,96} for Savonnier oolithic sandstone. The spikes at regular in- 
tervals for small L are not fluctuations, but rather a result of the anisotropy 
visible in Figure |TD| and Figure |TI|. 
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Figure 13: The family of two-dimensional LPDs ^2d{4''i L ') for L' G 
{1, 2, . . . , 240}. Note that the distributions are much broader than those 
for the overlapping sphere model shown in Figure |^. 
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tained using G\\ in (|^). The poor agreement for = 3 is to be expected, since 
G^ll 7^ G±, as seen in Figure 0. Hence, we cannot make an independent deter- 
mination of the three-dimensional variance from the two-dimensional data, 
and are unable to construct the mapping L (L) based on two-dimensional 
correlation functions. 

Nevertheless, a matching between two- and three-dimensional LPDs can 
be obtained based on A(L, L ) as discussed above. The results of the A- 
matching procedure are shown in Figure |1^ for a selection of (L, L )-pairs 
indicated in the figure. The agreement is of similar quality as that for the 
overlapping sphere model. 



7 Conclusion 

The present paper has investigated the dimensional dependence of local 
porosity distributions, which formalize the idea of length-scale dependent 
spatial averaging (coarse-graining). We have presented the first determina- 
tions of LPDs for real porous media. We have shown that in general, in 
the limit of large measurement cells(averaging regions), the two- and three- 
dimensional LPDs can be mapped onto each other by a rescaling of lengths 
L (L). This result was based on the central limit theorem and the concept 
of stationarity. We have found that the mapping L (L) extends below the 
asymptotic regime to much smaller length-scales. This result is important 
for both practical applications and because it agrees with the possibility of 
intermediate limit laws for LPDs based on generalizations of the central limit 
theorem 0. 

An equally important result of this paper is the finding that the map- 
ping L (L) can be determined from purely two-dimensional measurements. 
However, this result is restricted to homogeneous isotropic media. 

We illustrated our results with three examples. For the percolation model, 
the analytical form of the mapping L (L) between two- and three-dimensional 
LPDs is exactly known, and the agreement between mapped distributions 
becomes perfect. For the homogeneous isotropic overlapping sphere model, 
the agreement between mapped distributions is satisfactory for measurement 
cells of side L> ^ and improves with increasing L. The mapping L (L) can 
be successfully determined from two-dimensional observations alone for this 
medium. The three-dimensional pore-space reconstruction for the Savonnier 
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Figure 14: Two- and three-dimensional local porosity variances for the Savon- 
nier oolithic sandstone. The solid curves are variances determined directly 
from the distributions and the dashed curves are variances computed from 
(H) using the G\\ (r). 
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Figure 15: A(L, L')-matched two- and three-dimensional LPDs for the Savon- 
nier oohthic sandstone for a selection of (L, L )-pairs indicated in the figure. 
The two- and three-dimensional distributions are shown as dashed and solid 
curves, respectively. The spikes in Figure 12a are seen here in n^ai'P'i 8) 
(p = 1/8,2/8 and 3/8. 
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oolithic sandstone specimen is clearly anisotropic and shows large-scale het- 
erogeneities. The anisotropy of the reconstructed pore-space does not allow 
one to construct the mapping L (L) from two-dimensional observations. Nev- 
ertheless, the mapped curves from two- and three-dimensional measurements 
again show good agreement. 
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